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Abstract. The nature of the X(3872) enhancement is analysed in the framework of the Resonance- 
Spectrum Expansion, by studying it as a regular J PC = 1 ++ charmonium state, though strongly influenced 
and shifted by open-charm decay channels. The observed but Okubo-Zweig-Iizuka-forbidden p°J/ip and 
toj/ip channels are coupled as well, but effectively smeared out by using complex p° and to masses, in order 
to account for their physical widths, followed by a rigorous algebraic procedure to restore unitarity. A very 
delicate interplay between the D°D*°, p J/if), and uiJfip channels is observed. The data clearly suggest 
that the X(3872) is a very narrow axial- vector cc resonance, with a pole at or slightly below the D°D*° 
threshold. 



1 Introduction 



The A (3872) charmonium-like state was discovered in 2003 
by the Belle Collaboration [I], as a ir + ir~J/ip enhancement 
in the decay B ± — > K ± tt + it~ J/ip. The same structure was 
then observed, again in tt+tt" Jfip, by CDF II 0, DO [3], 
and BABAR g|. Moreover, CDF [5] showed that the tt+tt^ 
mass distribution favors decays via a p° resonance, imply- 
ing positive C-parity for the A(3872). The A(3872) has 
also been observed in the D°D°ir° and D*°D*° channels, 
by Belle [6] and BABAR [7j, respectively. More recently, 
CDF [5] measured the AT (3872) mass with even higher pre- 
cision, viz. 3871.61 ± 0.16 ± 0.19 MeV, with a width fixed 
at 1.34 ±0.64 MeV, while BABAR [9] presented evidence 
for the long-awaited uiJ/tp decay mode (also see Ref. [TUjh 
and a surprising preference for the 2 h assignment. We 
shall come back to the latter claim below. The A(3872) 
resonance is listed in the 2010 PDG tables [11] . with a 
mass of 3871.56 ±0.22 MeV, a width <2.3 MeV, and 1++ 
or 2 <~ quantum numbers |llU12j . 

On the theoretical side, the first to foresee a narrow 
1 ++ state close to the DD* threshold was Tornqvist [15] , 
arguing on the basis of strongly attractive one-pion ex- 
change for S'-wave meson-meson systems, which he called 
deusons. For further molecular descriptions and studies, 
see Ref. [14], as well as the reviews by Swanson [15] and 
Klcmpt & Zaitsev 16 . In Ref. [17], a few exotic model de- 
scriptions can be found, such as a hybrid or a tetraquark; 
also see the reviews [T51ll6j . For further reading, we rec- 
ommend the very instructive analyses by Bugg [TS] and 
Kalashnikova & Nefediev [T9l. 



Much more in the spirit of our own calculation is the 
coupled-channel analysis by Danilkin and Simonov |20) . 
which studies resonances and level shifts of conventional 
charmonium states due to the most important open and 
closed decay channels. We shall come back to their results 
below. 

As said above, according to the PDG [IT], the A(3872) 
is either a 1 ++ or a 2 h state, which implies 2 3 Pi or 1 1 Z?2, 
as other radial excitations would be much too far off (see 
e.g. Ref. [UJ). In the present paper, we only study the 1 ++ 
scenario, despite the conclusion by BABAR 9 , from the 
ojj/ifj mode, that 2 ^ is more likely. However, the latter 
assignment appears to be at odds with radiative-transition 
data [22 . For a further discussion of electromagnetic de- 
cays, see e.g. the molecular description of Ref. [25] , But 
more importantly, in all charmonium models we know of, 
the 1 1 Z?2 cc state lies well below 3.872 GeV, i.e., in the 
range 3.79-3.84 GeV (see e.g. Ref. [IS]). Our own bare 
1 1 D 2 state comes out at 3.79 GeV, just as the correspond- 
ing single-channel state in Ref. 20J . Now, the crucial point 
is that loops from closed meson-meson channels are always 
attractive [21]. Hence, since DD* at 3.872-3.880 GeV is 
the lowest Okubo-Zweig-Iizuka-allowed (OZIA) channel 
that couples to a 1 1 £>2 cc state, the coupled-channel mass 
shift will inexorably be further downwards (also see Ref. 
[25]). 

The goal of the present paper is to show that the mass 
and width of the V(3872), as well as the corresponding 
observed amplitudes in the D°D*°, p° J/ip, and ujJ/ip chan- 
nels, are compatible with a description in terms of a regu- 
lar 2 3 Pi charmonium state, though mass-shifted and uni- 
tarised via open and closed decay channels. 
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2 Resonance-Spectrum Expansion 

Sticking to the 1 ++ scenario, we employ again the Reso- 
nance-Spectrum Expansion (RSE) [25] , in order to couple 
one cc channel, with l c = 1, to several OZIA pseudoscalar- 
vector (PV) and vector-vector (VV) channels, just as in 
our recent preliminary study [27] of the A (3872). How- 
ever, we now also couple the OZI- forbidden (OZIF) p° J/tp 
and ujj/ip channels, to account for the bulk of the ob- 
served 7r + 7r - J/0 and ir + Tr~ir°J/ip decays, respectively. Al- 
though the former channel is isospin breaking as well, the 
extreme closeness of its central threshold at 3872.4 MeV 
to the A(3872) structure makes it absolutely nonnegligi- 
ble, despite a very small expected coupling. A complica- 
tion, though, is the large p width, which does not allow 
the p°J/tp channel to be described through a sharp thresh- 
old. Effects in the X(3872) from non-zero p and w widths 
were already estimated in Ref. [25]. We tackle this prob- 
lem by taking a complex mass for the p, from its pole 
position [21], and then apply a novel, empirical yet rig- 
orous, unitarisation procedure to the S-matrix, derived in 
Appendix [5] The analyticity and causality implications of 
complex masses in asymptotic states were already studied 
a long time ago [5U] . 

For consistency, we apply the same procedure to the 
to meson, despite the fact that its width is a factor 17.5 
smaller than that of the p. Nevertheless, the to width of 
about 8.5 MeV is very close to the energy difference be- 
tween the w J/0 and D°D*° thresholds, and therefore not 
negligible. Finally, we shall neglect the unknown small 
(< 2.1 MeV p]) D*° width [19], because of the relatively 
large error bars on the D°D*° data, though this width may 
have some influence on the precise A(3872) pole position. 
Nevertheless, reasonable estimates of the D*° width yield 
values clearly smaller than 100 keV 31 . so that its effect 
should be largely negligible as compared to that of the p 
and ui widths. 

Let us now proceed with our RSE calculation of a bare 
2 3 Pi (with n = 1, J = 1, L = 1, 5 = 1) cc state, cou- 
pled to a number of MM channels. The resulting closed- 
form T-matrix is given in Appendix [X] In Table [T] we 
list the considered PV and VV channels, including /r J/0 



Table 1. Included meson-meson channels, with thresholds and 
ground-state couplings. For simplicity, we omit the bars over 
the anti-charm mesons; also note that D* D* stands for the 
corresponding mass-averaged charged and uncharged channels. 



Channel 


G?(i c = l,n=0)) 


L 


Threshold (MeV) 


p°J/4> 


variable 





3872.406 -i 74.7 




variable 





3879.566 -i 4.25 


D°D*° 


1/54 





3871.81 


D°D*° 


5/216 


2 


3871.81 


D ± D* T 


1/54 





3879.84 


D ± D*^ 


5/216 


2 


3879.84 


D*D* 


5/36 


2 


4017.24 


DfD** 


1/54 





4080.77 


DfDt* 


5/216 


2 


4080.77 



and wJ/0. Besides the latter two OZIF channels and the 
also observed OZIA D°D*° channel, we furthermore ac- 
count for the OZIA PV and VV channels D ± D*^ 1 D* D* , 
and D S D* S , whose influence on the A(3872) pole position 
is not negligible, in spite of being closed channels. The 
D* S D* S channel, with threshold about 350 MeV above the 
A(3872) mass, we do not include. 

The relative couplings of the OZIA channels have been 
computed with the formalism of Ref. [32] . Couplings cal- 
culated in the latter scheme for ground-state mesons gen- 
erally coincide with the usual recouplings of spin, isospin, 
and orbital angular momentum. Moreover, for excited stat- 
es the formalism yields clear predictions as well, contrary 
to other appoaches. In Table [T] the squares of the ground- 
state (n = 0) couplings are given, which have to be multi- 
plied by (n + l)/4 ra for the S'-wave PV channels, and by 
(2n/5 + 1) /4™ for the others, so as to obtain the couplings 
in the RSE sum of Eq. ([3]). Also note that the two (closed) 
D*D* channels have been lumped together, with their av- 
erage threshold value and sum of squared couplings. As 
for the couplings of the p°J/ip and w J/0 channels, the for- 
malism of Ref. [32], based on OZIA decay via 3 Po quark- 
pair creation, cannot make any prediction. However, we 
know from experiment that the couplings of OZIF chan- 
nels are considerably smaller than those of OZIA channels. 
Moreover, isospin-breaking channels are even further sup- 
pressed. Thus, in the following we shall employ the values 
9p°j/iP = 0.07 x g D o D ,a and g uJ /^ = 0.21 x g D o D ,«, which 
correspond to effective relative strengths of 0.49% and 
4.41%, respectively, which seem reasonable to us. These 
values may also be compared to the corresponding rela- 
tive probabilites of about 0.65% (« 0.006/0.92) and 4.5% 
(« 0.041/0.92), respectively, employed in Ref. [23]. Fur- 
thermore, we shall also test coupling values twice as large, 
namely g p oj/^ = 0.14 x g D o D ,o and g uJ ^ = 0.42 x g D o D ,o. 
Note that our coupling for the isospin-breaking channel 
p J/0 is also in rough agreement with estimates from the 
rate of the observed [TT] isospin- violating oj — » tt + tt~ de- 
cay, which amounts to about 1.5% of the total width. An- 
other difference between OZIA and OZIF channels is the 
average distance (see Eqs. ( Tj2 ) ) at which a light qq 
pair is created before decay, which in the OZIA case we 
believe to take place in the core region and in the OZIF 
case more in the periphery. Thus, we employ a larger 
value for r\ = r p aj/^, — r u M than for ro, the single 
radius used for all OZIA channels. Concretely, we take 
r = 2 GcV- 1 ~ 0.4 fm and n = 3 GcV" 1 ~ 0.6 fm, 
while we also test the case r-y = r$. 

For the bare cc energy levels En c ^ in the RSE sum of 
Eq. we take an equidistant harmonic oscillator (HO), 
as in previous work (see e.g. Ref. [33]). The corresponding 
charm quark mass and HO frequency are once again kept 
fixed at the values m c = 1562 MeV and u = 190 MeV. 
The only parameter we adjust freely is the overall coupling 
constant A in Eqs. ( 1|2 1, which is tuned to move the bare 
2 3 Pi state from 3979 MeV down to the D°D*° threshold, 
requiring a A value of the order of 3, i.e., not far from the 
values used in e.g. Refs. [551134"] . At the same time, the bare 
1 3 Pi state shifts from 3599 MeV down to about 3.55 GeV, 
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though depending quite sensitively on the precise form of 
the used subthreshold suppression of closed channels [33] . 
Anyhow, for the purpose of the present study, an accurate 
reproduction of the \ c i(lP) mass of 3511 MeV is not very 
relevant. 



3 X(3872) poles and amplitudes vs. data 

In Table [2l we give some pole positions in the vicinity of 
the D°D™ and p°J/ip thresholds, with the chosen values 



Table 2. Pole positions of the dots and stars in Fig. [l] In all 
cases, n = 3.0 GeV -1 . Note that the OZIF couplings g p oj/^ 
and m are given relative to the coupling of the OZIA D°D*° 
channel. 



Label 


A 


9 P °j/^ 




Pole (MeV) 


1 


3.028 


0.07 


0.21 


3872.30 -i 0.71 


2 


3.066 


0.07 


0.21 


3871.83 -i 0.40 


3 


3.083 


0.07 


0.21 


3871.56 -i 0.11 


a 


2.981 


0.14 


0.42 


3872.30 -i 0.75 


b 


3.017 


0.14 


0.42 


3871.82 -i 0.48 


c 


3.033 


0.14 


0.42 


3871.57 -i 0.28 



of A and r\ . In Fig. [T] third-sheet pole trajectories in the 
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Fig. 1. Pole trajectories for ri=3.0 GeV (solid curves) and 
ri=2.0 GeV -1 (dotted curves); g p aj^/g D o D ,a = 0.07 and 
9wJN>/9d°d*o = 0.21 (upper two curves), g p0 J/0 / 'g D o D ,a = 0.14 
and g^j/^/goOD* = 0.42 (lower two curves). Note that the CM 
energy E is relative to the D°D*° threshold in all figures. Also 
see Tabled 



complex energy plane (relative to the D°D*° threshold) 
are plotted, as a function of A, with the pole positions 



of Table [2] marked by bullets and stars. The solid curves 
represent the case rj = 3.0 GeV -1 , while the dotted ones 
stand for r\ = 2.0 GeV -1 , showing little sensitivity to 
the precise decay radius. Figure [l] shows that the X(3872) 
resonance pole may come out below the D°D*° threshold 
with a nonvanishing width, which is moreover of the right 
order of magnitude, viz. < 1 MeV. The recent CDF 8 
mass determination of the X(3872) might suggest that 
the pole positions '3' or 'c' (see Table [2] and Fig. [I]) are 
favored. However, one should realise that the differences 
amount to mere fractions of an MeV, while experimental 
uncertainties are at least of the same order. 

Now we compare the corresponding D°D*° amplitudes 
to Belle [3S] data , for the six cases labeled T, 2, 3', and 
'a, b, c' in Tableland Fig.[l] The results are depicted in 
Figs. [2] and [3j respectively. Note that we allow for an arbi- 



40- 



30- 



20- 



10- 




2 4 6 

E (MeV) 

Fig. 2. D°D*° elastic amplitude for poles 1, 2, 3 in Tableland 
Fig. [l] arbitrarily normalised data are from Ref. [35]. Elastic 
T-matrix elements follow from Eqs. (J1J-J3J) ; k is on-shell relative 
momentum. 



trary normalisation of the data, which is inevitable as we 
are dealing with production data, which cannot be directly 
compared with our scattering amplitudes, also because of 
the finite experimental mass bins. From these figures we 
see that the best agreement with data is obtained in case 
'2', though 5 out of the 6 curves pass through all error 
bars. Nevertheless, in view of the large errors, one should 
be very cautious in drawing definite conclusions on the 
precise pole position as well as the preferred OZIF cou- 
plings g p oj/^ and g^j/ip. 

Next we show, in Fig. [1J the elastic amplitudes in the 
p J/ip and ujJ/ip channels, corresponding to the pole po- 
sitions 1, 2, 3, i.e., for the smaller values of the OZIF 
couplings. We see that both amplitudes are very sensitive 
to the precise pole position, which is logical, as the OZIF 
channels couple much more weakly to cc than D°D*° 1 so 
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Fig. 3. As in Fig. [2J but now for poles a, b, c. 
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Fig. 4. p a J/ip (left) and wJ/ty> (right) elastic amplitudes for 
poles f, 2, 3. Also see Fig. [l] and Table(2] 



that the latter channel will strongly deplete the former 
ones, as soon as it acquires some phase space. This is 
in line with our analysis in e.g. Ref. |36j . Also note the 
strongly cusp-like structure of the amplitude in the cases 
2 and 3 for p° J/ip, and 3 for a; J/0, which is a manifestation 
of the depletion due to the opening of the D°D*° channel. 
Such a cusp makes the experimental determination of the 
Jf(3872) width very difficult. 

In Fig. [5] we take a closer look at the u) J/ip and p°J/ip 
amplitudes, in particular how they compare to one an- 
other. Now, the effective strength of the uJ/ip elastic T- 
matrix element is 9 times that of p°J/ip, as its coupling 
has been chosen 3 times as large (see Tableland Eqs. ([!]- 
[3|). For the corresponding square amplitudes plotted in 
Fig. [5j this amounts to a factor as large as 81. However, 
the central cjJ/0 threshold lies more than 7 MeV above 
that of p°J/ip, while the full w width is only 8.49 MeV. 
On the other hand, the central p°J/ip threshold lies much 
closer to D°D*°, while the large physical p width strongly 
boosts the associated amplitude, as demonstrated below. 
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Fig. 5. p°J/ip and u)J/ip elastic amplitudes, for poles 2 (left) 
and 3 (right). Also see Fig. [I] and Table [5] 



Qualitative arguments in agreement with our calculation 
were already presented in Ref. [37] . These effects make the 
maximum uiJ/ip square amplitude to be only a factor 3.5— 
4 larger than that of p°J/tp, both in case 2 and 3, as can 
be read off from Fig. [5] Moreover, at the precise energy 
of the respective pole position, the two amplitudes are al- 
most equal in size. Therefore, the observed branching ratio 
£(AT(3872) uiJ/ip)/B(X(3872) -> n+n-J/ip) ~ 1 0HU] 
is compatible with the present model calculation. 

Finally, in order to study the effect of using a com- 
plex mass for the p° in the p°J/ip channel, we vary the p 
width from 0% to 100% of its PDG QI] value and plot 
the corresponding amplitudes in Fig. [6j We see that the 
maximum \T\ 2 increases by almost 3 orders of magnitude 
when going from the 0% case (dotted curve in left-hand 
plot) to the 100% case (solid curve in right-hand plot). 
Furthermore, the 0% curve only starts out at the central 
p°J/ip threshold, of course. Thus, it becomes clear that 
no realistic description of the p°J/ip channel is possible 
without smearing out somehow its threshold, so that its 
influence kicks in before the D°D*° channel opens and 
depletes the signal. Naturally, similar conclusions apply 
in principle to the uJ/ip channel, though there the effects 
are less pronounced because of the small uj width and the 
somewhat higher threshold. These results show that our 
unitarisation procedure for complex masses in the asymp- 
totic states performs as expected in accounting for thresh- 
olds involving resonances. 

To conclude our discussion, we should mention that 
our results are qualitatively in agreement with those of 
Danilkin & Simonov [5D], in the sense that a single res- 
onance pole originating from the 2 3 Pi cc state is capa- 
ble of describing the AT (3872) data. However, we disagree 
with their conclusions on the 2 3 Po state. In an earlier, 
single-channel description |38j. we found a resonance at 
3946 MeV with a width of 58 MeV, and we do not be- 
lieve a detailed multichannel calculation will change these 
values dramatically. Thus, the listed A"(3945) [11] reso- 
nance, with mass 3916 MeV, width 40 MeV, and positive 
C-parity, appears to be a good candidate. As for the 2 1 P\ 
state, the AT(3940) QI] resonance, with mass 3942 MeV, 
width 37 MeV, and principal decay mode DD* , seems the 
obvious choice. With the old Z(3930) meanwhile identified 
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Fig. 6. p°J/ip elastic amplitude for reduced p° width. Left: 0% 
(dots), 1% (dashes), 5% (full); right: 10% (dots), 50% (dashes), 
100% (full). Studied case: pole 2 in Table [2] Also see Fig.[l] 

as the 2 3 P2 (Xc2(2P) [H]) state, we might so understand 
all 4 charmonium states in the range 3.87-3.95 MeV. 



4 Summary and conclusions 

Summarising, we have investigated the 1 ++ charmonium 
scenario for the AT (3872) resonance, by analyzing in de- 
tail the influence of the D°D*°, p°J/ip, and coJ/tp channels 
on pole positions and amplitudes. In order to describe the 
latter OZIF channels in a realistic way, we have used com- 
plex masses for the p° and oj, and then restored unitarity 
of the ^-matrix by a new and rigorous algebraic procedure, 
albeit physically heuristic. It is true that the redefined S- 
matrix may have some unusual analyticity properties |30] , 
but in our amplitudes no sign was found of any nearby spu- 
rious singularities. Moreover, the behaviour of the p°J/ip 
amplitude as a function of the p° width gives us confi- 
dence in our approach. Concretely, we have shown that 
our scenario is compatible with the D°D*° and ir + ir~ J/i/> 
data, with a single resonance pole on top of or slightly 
below the D°D*° threshold. Moreover, our treatment of 
the p° J/0 and uiJ/i/j channels has proven compatible with 
the observed branching ratio of these decays. 

Thus, the data do not seem to require a molecular or 
tetraquark interpretation of the AT(3872), also in view of 
so far unobserved charged partner states. Neverthe- 
less, only further improved measurements and theoretical 
calculations will in the end allow to draw a definitive con- 
clusion on the scenario preferred by nature. 

In conclusion, we must stress that the V(3872), what- 
ever its assignment, is an extraordinary structure, because 
of its coincidence — to an accuracy of less than 1 MeV — 
with the central thresholds of the principal decay modes. 
This circumstance is at the same time a blessing and a 
curse. To start with the latter, no model can ambition 
to quantitatively describe the V(3872) with present-day 
state-of-the-art in strong interactions, while experiment 
will have an extremely hard time to reduce the bin sizes 
to less than 1 MeV and simultaneously keep statistics suf- 
ficiently high. On the other hand, with a strengthened 
effort of both theory and experiment, a wealth of knowl- 
edge on charmonium spectroscopy and strong decay - 



OZI-allowed as well as OZI-forbidden — may be gathered 
by further studying this fascinating resonance. 

We are indebted to Prof. D. V. Bugg and J. Segovia for very 
useful comments. This work was supported by the Fundacao 
para a Ciencia e a Tecnologia of the Mmisterio da Ciencia, 
Tecnologia e Ensmo Superior of Portugal, under contracts nos. 
CERN/FP/83502/2008 and CERN/FP 109307/2009. 



A Multichannel T-matrix 

The off-energy-shell RSE T-matrix reads [26 1 133 ] 

T^-^ip^-E) = -2A 2 7^^jl i ( K r l ) x 

N 



\i0 fj 



(1) 



with the diagonal loop function 

fyj(kj) = 2/A-// ,/.-,r, f L . {kjrj) h^'/ (kjTj) S l3 , (2) 
and the RSE propagator 



n, 



E 



p(ic) 



(3) 



Here, the RSE propagator contains an infinite tower of 
s-channel bare qq states, corresponding to the spectrum 
of an, in principle, arbitrary confining potential. Further- 
more, A is an overall coupling, is the decay radius of 

meson-meson (MM) channel i, En"^ is the discrete en- 
ergy of the n-th recurrence in the qq channel with angu- 
lar momentum l c , g 1 ^ ^ is the corresponding coupling to 
the i-th MM channel, //, the reduced mass for this chan- 
nel, Pi the off-shell relative momentum , and j l L .{Pi) an< ^ 

(kjVj) the spherical Bessel and Hankel functions of 
the first kind, respectively. Note that /ij, pi, and fej (the 
on-energy-shell relative momentum) are defined relativis- 
tically (see e.g. Ref. (ID], Eqs. (16,17)). Also notice that, 
for a J PC = 1 ++ cc system, there is only one confined 
channel, viz. with l c — 1. 

Bound states and resonances can be found by search- 
ing for zeros in the determinant of the matrix ( 1 — Q TV) 
in Eqs. (Uj[3]), employing a complex Newton's method. On 
the other hand, elastic and inelastic amplitudes are ob- 
tained by taking the on-shell values of the corresponding 
matrix elements of T. 



B Redefining the 5-matrix 

It is straightforward to show that the S'-matrix S(E) 
[1 + 2iT(E), where T(E) is the on-energy-shell restriction 
of the multichannel T-matrix in Eqs. (JTJ-J3]) , is unitary and 
symmetric, when limited to open channels and real en- 
ergies. However, it is also easy to see that, for complex 
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masses and so complex relative momenta, the unitarity of 
5 is lost, though not its symmetry. The latter property 
can be used to redefine the physical 5-matrix. 

Since 5 is always a symmetrix matrix, it can be de- 
composed, via Takagi |3Tj factorisation, as 

5 = VDV T , (4) 

where V is unitary and D is a real nonnegative diagonal 
matrix. Then we get 

5 f 5 = (V T )^DV^VDV T = (V T ) f D 2 V T = U^D 2 U, (5) 

where we have defined U = V T , obviously unitary, too. So 
the diagonal elements of D = y/US^SW are the square 
roots of the eigenvalues of the positive Hermitian matrix 
5^ S, which are all real and nonnegative. Moreover, since 
5=1 + 2iT is manifestly nonsingular, the eigenvalues of 
5^5 are even all nonzero and U is unique. Thus, we may 
define 

5' ee SrfD-^U . (6) 
Then, using Eq. Q and V = U T , we have 

5' = XJ T DXJU^D- Y U = U T U , (7) 

which is obviously symmetric and, as 

(U T Uy = t/t([/t)T = U-lp-lf = {U T ur l f (8) 

also unitary. So 5' has the required properties to be de- 
fined as the 5-matrix for a scattering process with complex 
masses in the asymptotic states. 
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